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The valence-band spectrum of the ferromagnetic nickel is calculated using the LDA+DMFT 
method. The auxiliary impurity model emerging in the course of the calculations is discretized 
and solved with the exact diagonalization, or, more precisely, with the Lanczos method. Particular 
emphasis is given to spin dependence of the valence-band satellite that is observed around 6 eV 
below the Fermi level. The calculated satellite is strongly spin polarized in accord with experimental 
findings. 
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I. INTRODUCTION 

The electronic structure of transition metals has been 
intensively studied for a number of decades. Notwith- 
standing, certain aspects of the electron behavior in these 
materials are still not completely understood. Comparison 
of experimental findings with effective one-electron band 
theories have indicated that a more thorough treatment 
of quantum many-body effects is necessary to accurately 
describe the physical reality. 

A prototypical metal displaying pronounced electron 
correlations is the ferromagnetic nickel, where the one- 
particle spectrum obtained using the local-density approx- 
imation (LDA) to the density-functional theory (DFT) 
noticeably departs from the spectra measured in photoe- 
mission experiments. The calculated 3d bandwidth as 
well as the exchange splitting are too large.^^ Moreover, 
the LDA completely misses the satellite feature located at 
approximately 6 eV below the Fermi level. This satel- 
lite was originally attributed to plasmon excitations,'^ but 
an alternative view was soon proposed,!^ according to 
which the satellite is a result of a correlated state of two 
Sd holes localized in a single atom. The latter picture 
is supported by the observed resonant enhancement of 
the satellite, during which a second scattering channel 
involving 3p electrons and ending in the same two-hole 
final state becomes active.l^Win! The correspondence be- 
tween the satellite and the localized two-hole states can 
be explicitly visualized in simplified finite- sized models 
that allow for an exact many-body solution.EnH 

A more quantitative description of the electron correla- 
tions in nickel can be achieved by incorporating a selfen- 
ergy into the LDA or Hartree-Fock bandstructure.'^^'^ 
Usually, the selfenergy is assumed local, that is, wave- 
vector independent. The most sophisticated local self- 
energy isprovided by the dynamical-mean-field theory 
(DMFT}i3 that maps the problem of interacting lattice 
electrons onto an impurity model where the interactions 



are restricted to a single lattice site. The combination 
of LDA and DMFT (the so-called LDA-hDMFT method) 
was applied to the electronic structure of nickel several 
times in the past, using different methods to solve the 
auxiliary impurity model.'^^l^ A reasonable description 
was achieved employing the Hirsch-Fye quantum Monte 
Carlo (QMC) method as the impurity solver.C^The QMC 
methods have many merits. In particular, they are consis- 
tently accurate regardless of the strength of correlations 
in the system. But they have weaknesses too. The QMC 
calculations of the one-particle spectral functions involve 
a numerical continuation from the imaginary time to the 
real frequencies, a procedure with a limited resolution 
especially at higher binding energies. Additionally, the 
QMC algorithm used in Ref. [17] is limited to a diagonal- 
only Coulomb interaction. This truncation breaks a subset 
of symmetries characterizing the full Coulomb operator, 
which can lead to undesirable side effects. 

In this paper we solve the auxiliary impurity model of 
the LDA+DMFT by means of the Lanczos method. This 
strategy involves a discretization of the impurity model 
which represents an obvious limitation on the achievable 
accuracy. The sources of errors in this approach are, 
however, very different from those in the QMC method 
and the two impurity solvers can thus offer complementary 
information. Using the Lanczos method, the one-particle 
Green's function can be evaluated directly anywhere in the 
complex plane without resorting to any extrapolation. It 
is also straightforward to compare the full and truncated 
Coulomb operators, and we make this comparison in 
Secini 

II. METHOD 

We start from the bare electronic structure of Ni ex- 
pressed in terms of a tight-binding LMTO modeP^l con- 
taining 4s, 3d and Ap electronic states. The one-electron 
Hamiltonian H{k) is obtained as a solution of the local- 
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density approximation and the correlations beyond this 
approximation are accounted for by a local selfenergy E 
acting in the subspace of the d orbitals. The selfenergy 
is spin polarized whereas the Hamiltonian H(k) is spin 
independent. Taking H(k) from the spin-polarized LDA 
is also possibfejalthough this route was found as less 
accurate earlierpSJ 

The selfenergy S is constructed with the aid of an im- 
purity model defined by a Hamiltonian ^imp = -^^hnp + ^ 
that describes a single d shell hybridized with a sea of aux- 
iliary conduction electrons. These auxiliary electrons, of- 
ten referred to as the bath, model the environment around 
the d shell in the actual nickel lattice. The Coulomb 
interaction U acts only among the d orbitals and the 
one-particle part H^^^ has the form 
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), (1) 



where d\^^ creates an electron in the d shell and c^^^ 
creates a conduction electron in the bath. The index m 
runs over Cg = {x^ — y^, z^} and t2g = {xy, xz, yz} states, 
and a E {t, i} labels spin projections. The hybridization 
parameters Vkma- couple only those impurity and bath 
levels that carry the same indices m and tr, and hence 
the cubic symmetry and the electron spins are preserved. 

Provided we can solve the interacting impurity model, 
the sought for selfenergy E is obtained as 



^ ~ '^imp [^imp] ^imp [-^imp] 



(2a) 



where Gi^p[H] represents the Green's function matrix in 
the d-orbital subspace evaluated for a general impurity 
Hamiltonian H. The matrix Gimp [^/mp] i which we will 

denote as G for short, is usually referred to as the bath 
Green's function. 

So far, we have not specified how the parameters en- 
tering the Hamiltonian TJimp should be determined. The 
missing link to the original lattice electrons is provided 
by a condition 



Gimp[^imp] =G[i?(k),E] 



(2b) 



that equates Gimp to the local d-orbital Green's func- 
tion G evaluated in the lattice. The right-hand side of 
Eq. ( 2b I can be expressed as a momentum sum over the 



first Brillouin zone 



(3) 



where / stands for the identity operator and the chemical 
potential is chosen such that the As-3d-4p space holds 
ten electrons per Ni atom. 



Equations (|2| define the dynamical-mean-field approxi- 
mation. They are iteratively solved for S and H^^^ while 

iJ (k) and U act as inputs. The most involved part of this 
procedure is the solution of the multi-orbital impurity 
model. A number of approximations of varied accuracy 
have been used to find this solution in the context of 
the DMFT. Here we discretize the impurity Hamiltonian 
-ffimp and then solve the resulting finite-sized cluster He 
essentially exactly by means of the Lanczos method. This 
strategy was suc cessfu lly applied to the DMFT equations 
for the repulsiv^ISESl ^nd attractiv^^i! single-band Hub- 
bard models as well as for realistic multi-band problems.'^ 
The discretization i/imp — > He amounts to a replacement 
of the infinite sums (integrals) over k in Eq. ([T]) with short 
finite sums. In our particular case, the index k takes only 
two values, that is, each impurity orbital is connected to 
just two bath orbitals. 

The parameters of the discretized Hamiltonian {ema, 
^kma and Vkma) are expressed as functions of H(k) and 
E with the aid of the relation 



g-' EE Grm'p [ff (°)] « G-1 [H{k), ±]+±^g- 



(4) 



which is just a rearranged form of Eqs. ([2]). At this point 
it is necessary to specify in what sense the discrete bath 
Green's function Qc approximates the continuous function 
g , that is, what is the precise meaning of the symbol ~ 
in Eq. Q. It has become customary to minimize some 
distance between Gdz) and g{z) defined on the Matsubara 
frequencies z = kon- A p articularly convenient choice is a 
least-squares fit |^^ l ^'^* ^^for instance 



min y 
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for each m and a. For reasons that will be discussed 
later, we do not follow this fitting route but adopt an 
alternative approach instead. We obtain the parameters 
of the discretized Hamiltonian from the requirement of 
coincidence of the high-frequency asymptotics of Qc{z) 
and G{z)^^ 

The cluster Green's function Qc{z) can be written in 
an explicit forrrP^I 



E 



(6) 



whose expansion in powers of 1/z reads as 
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The continuous Green's function Qiz) can be expressed in 
terms of the density of states g(z), and the coefficients of 
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the expansion in powers of 1 j z are then given as moments 
of this density of states, 



(8) 



Mn-l 



With two bath orbitals per each impurity d orbital we 
have five parameters in iJc"^ that carry the same indices 
m and cr, and thus we can match Eqs. ([7| and ^ up to 

We are mostly interested in the ground-state properties 
and in the one-particle spectrum in, say, the first 10 eV 
below the Fermi level. In the course of our calculations 
we observed that the discretization procedure defined 
by Eqs. (|4])-([8| often placed some of the bath energies 
Cfemo- quite high above the Fermi level far outside the 
energy window of our interest. That by itself would not 
be an issue if it did not lead to an unphysical stabilization 
of a non-magnetic solution. (See Appendix |b] for an 
illustration and further discussion of the effect). In order 
to suppress this undesirable behavior, we modify the 
definition of the moments M„ to 



M„ = 



gma{e) de 



(9) 



The lower cutoff is a purely technical matter; it is set to 
ci = —9 eV, that is, below the 4s band. The upper cutoff 
avoids the unphysical solution by not allowing the bath 
orbitals to drift to high energies. The results presented 
in Sec. |III| were obtained with = 2 eV. The possibility 
to straightforwardly prevent the non-magnetic state with 
the aid of the upper cutoff is the main reason why 
we opted for the bath discretization by means of the 1/z 
expansion instead of the more frequently employed fitting 
on the Matsubara axis. We have not succeeded in finding 
a suitable modification of the fitting function, Eq. ([s]), 
that would reliably eliminate the non-magnetic solution. 

The last component of the cluster Hamiltonian He is the 
Coulomb interaction in the d shell. We use the spherically 
symmetric form 



U 
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where the matrix Umm'm"m"' is parametrized by the 
Slater integrals Fq = 2 eV, F2 = 8.2 eV and = 5.2 eV. 
These numerical values correspond to Coulomb [/ = 2 eV 
and exchange J = 0.95 eV. The contribution to Eq. (10) 
proportional to J7h represents a rigid shift of the impu- 
rity levels downward, Cma — >■ ^ma — C/h, and accounts for 
the fact that the d-d Coulomb interactions are already 
partially included in the LDA Hamiltonian H(k) in the 



form of a static mean field. Several formulas have been 
proposed to express the Hartree-like double-counting po- 
tential U}i in terms of the occupation numbers of the d 
orbitalSjI^SHsSl but we treat Un as a free parameter sim- 
ilarly to Ref. |3TJ since neither of the standard choices 
leads to satisfactory results. 

The need for "undressing" the LDA quasiparticles is 
one of the reasons why we prefer to build the many- 
body description on the top of the spin-restricted LDA 
bandstructure. If we started from polarized bands, the 
Hartree potential Un would be polarized too, which would 
introduce an extra complexity to the problem. The double 
counting would have to be spin dependent also in the 
LDA-I-DMFT implementations that take into account the 
feedback of the selfenergy on H(k)W^ 

Once the cluster Hamiltonian is fully specified, the 
one-particle Green's function Gc = Gimpi-ffc] for individ- 
ual d orbitals can be calculated. We employ the band 
Lanczos method'^'' that allows for a simultaneous evalu- 
ation of all relevant matrix elements at once. Off-diagonal 
elements are directly accessible too, although this func- 
tionality is not used in the application at hand. For the 
purpose of the Lanczos method, Gc is decomposed in the 
following forrrPH 

G'„.aiz)^^[G^.aiz)+G<,{z)], (11) 

where the two parts are 

a 
a 

The sums over the many-body eigenstates \a), 
Hc\a) = Ea\a), represent grandcanonical aver- 
ages with the chemical potential equal zero, and 
Z = J2a Q^^^" stands for the corresponding partition 
function. The calculations are performed at low temper- 
ature UbT = 1/(3 — 0.01 eV so that only the ground 
state contributes to the sum over a most of the time. The 
eigenstate-eigenvalue pairs including all their degeneracies 
are found using the implicitly restarted Lanczos method 
as implemented in the ARPACK software package.!^ 

III. RESULTS AND DISCUSSION 
A. Ground-state properties 

First we examine selected characteristics of the ground 
state and use them to estimate the double-counting po- 
tential [/h- Figure [T] shows the number of electrons in 
the d orbitals — rid-]- + ridi and the spin polarization of 
these orbitals rrid — rid^ — iidi- The data calculated in the 
lattice and in the discretized impurity model are plotted 
side by side. They differ despite the DMFT iterations 
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FIG. 1. (color online) The occupation of the d orbitals 
Ud (top) and the spin polarization md (bottom) plotted as 
functions of the double-counting potential Uh- Empty symbols 
correspond to the cluster Green's function Gc, full symbols to 
the lattice Green's function G. 



being converged in the sense that the cluster Hamilto- 
nian He no longer changed in the successive steps. The 
differences would vanish if we solved the full continuous 



impurity model, since Eq. ( 2b ) would be exactly fulfilled 
in that case. 

It turns out that rid and rn^ depend only weakly on the 
double-counting potential [/jj when the latter is larger 
than approximately 14.5 eV. Below 14.5 eV the trend 
changes and the cluster quantities depart substantially 
from their lattice counterparts. Based on this observation 
we consider Un below 14.5 eV as inappropriate. We note 
in passing that the double counting in the so-called fully 
localized limit™ c/^fll) ^ ^^^^ _ ^^^^ _ j^^^ _ ^^^^ 

equals 13.2 eV for n^j = 9 and it is thus more than 1 
eV too small to be applicable in our case. The so-called 
around mean- field forrrP^ of C/h provides an even smaller 
value. 

The experimentally determined magnetization of the 
fee nickel is approximately 0.6 fj,^ per atom.l^S! Our calcu- 
lations slightly underestimate this quantity even though 
the cluster solution, from which the spin-dependent self- 
energy is extracted, displays the maximal polarization 
characterized by md = 5 — ti^^. 

The number of d electrons cannot be unambiguously 
defined in a solid and as such it does not represent a 
particularly useful measure of quality of our ground state. 
The d-band filling in nickel is often estimated as 9.4 per 
atom based on the measured magnetic moment and the 
assumption of the maximal d-shell polarization,!^ but 
reliability of this estimate is limited. 



tuWU 




■e 

o 



full U \ 






ll D I 

!i — bu 


FN 








■e 

o 



FIG. 2. (color online) Spin-resolved d-orbital spectral function 
of bulk Ni obtained with f/n = 15 eV. The middle panel 
corresponds to the full Coulomb vertex, the bottom panel 
to the truncated vertex. The atomic d? transitions are 

displayed at an arbitrary scale and position in the middle panel 
for comparison with the shape of the satellite. The top panel 
shows the d-orbital occupation in the first 5000 many-body 
final states corresponding to the cluster Hamiltonian with the 
full vertex. 



B. Valence-band spectrum 

We find that the one-particle spectra corresponding 
to the double-counting potential [/h in the range 15.0 ± 
0.5 eV are only barely distinguishable. Figure [2] shows 
the d-orbital spectral function IniX^m ^ iO)]/^ 

for [/h = 15 eV. The calculated spectrum is relatively 
disappointing: the width of the main band (w 4 eV) as 
well as the exchange splitting are nearly identical to those 
obtained with the spin-polarized LDA, and thus share the 
same poor agreement with experiments. The symmetry- 
resolved exchange splitting at the Fermi level is given 
directly by the selfenergy and reads as 

Ee,T(^F)-Se^i(£^F)«0.3eV, (12a) 
E,,,^(i?F) - ^t^^^E^) « 0.8 eV. (12b) 

The d states near the Fermi level have predominantly 
the t2g character, which results in the apparent exchange 
splitting of 0.6 eV that is visible as a shift between the 
top of the valence bands for the minority and majority 
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spins (indicated with arrows in Fig. [2]) . Appendix |A] 
indicates that the band width and the exchange sphtting 
could possibly be improved if the impurity model was 
discretized using Eg. ([sj) instead of Eqs. ^ and ([s]), but, 
as discussed in Sec. that approach is not sufficiently 
reliable. 

We identify the spectral features below —4.5 eV as the 
"6 eV satellite". It is strongly spin polarized in agree- 
ment with spin-resolved photoemission experiments."^ In 
our calculations, the energy-integrated spectral weight 
is about three times larger for the majority spins than 
for the minority spins. Furthermore, the minority-spin 
states are located at reduced binding energies, which was 
also observed experimentally.ES The calculated character- 
istics of the satellite corroborate its explanation based on 
transitions from the spin-polarized cP initial state to the 
final states. An illustration of such atomic spectral 
lines is added to Fig. [2] for comparison. The singlet final 
states ^D, and exhibit a complete majority-spin 
polarization and lie deeper, the triplet states "^F and '^P 
carry a partial polarization in the opposite direction and 
lie shallower. 

This simplified description of the satellite should not 
be taken too literally, however, at least not within our 
computational scheme. We have calculated the d-orbital 
occupation corresponding to the final states in our 
discretized impurity model, the results are aligned with 
the lattice spectral function in Fig.|2] Although Ud indeed 
decreases as the binding energy increases, it is still con- 
siderably larger than eight in the satellite region where 
contributions from states with > 8.5 are not an excep- 
tion. This enhancement of is due to impurity-bath 
hybridization as discussed at the end of Appendix [Bj It 
is possible that Ud is somewhat overestimated as a result 
of the sparse discretization of the bath. 

As mentioned earlier, our calculations are rather insen- 
sitive to a choice of the potential J7h as long as it exceeds 
a threshold of approximately 14.5 eV. For smaller Uh the 
impurity orbitals in the cluster start to depopulate, which 
is accompanied by an increased intensity of the satellite. 
This result is in accord wit h exp eriments on alloys of Ni 
with electropositive metals.BSHSl 

Finally, we compare spectral functions calculated with 
two versions of the Coulomb operator: the full spherically 
symmetric vertex discussed so far, and the diagonal-only 
vertex employed in the Hirsch-Fye QMC method.*^ Fig- 
ure [2] shows that the simplification of the interaction 
has virtually no effect on the main d bands. The satel- 
lite, on the other hand, is visibly modified. The tail of 
the majority-spin spectrum does not extend as deep as 
with the full interaction, and the minority-spin satellite 
is shifted to smaller binding energies. It is possible that 
this slight shift in conjunction with the low resolution 
of the maximum-entropy method leads to a merger of 
the satellite with the main band for the minority spins, 
resulting in a fully spin-polarized satellite reported in 
Ref.HZl 



IV. CONCLUSIONS 

We have investigated the valence-band spectra of the 
ferromagnetic nickel within the LDA-(-DMFT framework. 
The auxiliary impurity model was discretized and then 
solved using the Lanczos method. The valence-band 
satellite and its spin dependence was reproduced in good 
agreement with spin-resolved photoemission experiments. 
The many-body renormalization of the Sd bands as well 
as the exchange splitting are found to be sensitive to the 
details of the bath discretization, which indicates that ten 
orbitals are probably not enough to approximate the bath 
to a satisfactory accuracy. The diagonalization method as 
employed in this paper is adequate for recovering features 
of atomic origin located at high binding energies but it is 
apparently too crude to capture the expected modification 
of the Fermi-liquid parameters at low binding energies. 
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Appendix A: Comments on the bath discretization 

In this appendix we return to the bath discretization 
using the least-squares fit of the bath Green's function 
with the functional form of Eq. Q at the Matsubara 
frequencies. In the case of five orbitals in the bath, the 
self-consistent LDA+DMFT calculations are very stable 
and the resulting spectral function is shown in the Fig. [3] 
Although the satellite at high binding energies is not de- 
scribed as well as previously (Fig. [2]), the renormalization 
of the main valence band and the exchange splitting come 
out as more reasonable. The width of the valence band is 
approximately 3 eV and the exchange spitting is about 
0.3 eV, both of which are close to the photoemission 
experiments.'^ The fitting at the Matsubara frequencies 
is more sensitive to the behavior near the Fermi level than 
the method of moments, Eqs. Q and ([s]), and hence it 
gives a finer control over the low-energy spectral features. 
Unfortunately, due to the effects discussed in the next 
appendix, it is not easy to converge the calculations to 
the correct state for ten orbitals in the bath. Neverthe- 
less, we believe that it will eventually be possible to find 
an optimal way of fitting the bath Green's function that 
would lead to a reasonable description of both low- and 
high-energy parts of the spectrum simultaneously. 
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FIG. 3. (color online) Spin-resolved d-orbital spectral function 
of bulk Ni that was obtained with the bath discretized using 
the least-squares fit at the Matsubara frequencies (five orbitals 
in the bath). 




-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 

FIG. 4. (color online) Spin polarization of the d shell (top) 
and the total number of electrons in the whole cluster ntot 
(bottom) as functions of the bath-level position eb. 



Appendix B: Magnetism in the finite cluster 

Here we take a closer look at the issues that necessitated 
the introduction of the cutoffs in Eq. (|9]). To illustrate 
the problem, we use a simpler impurity model than that 
employed in Sec. |III| — we reduce the cluster to contain 
only one bath orbital per each impurity orbital and we 
assume spherical instead of cubic symmetry. Furthermore, 
we implement the Hamiltonian parameters used in Ref . I12[ 
which gives us the opportunity to relate our calculations 
to this earlier study of electron correlations in nickel. The 
Slater integrals are Fo = 3.5 eV, F2 = 9.79 eV and F4 = 
6.08 eV, and the impurity-bath hopping is Vkma = 0.7 
eV. The bath-level position tkma = is treated as a 
free parameter and the double-counting potential C/h is 
determined such that there are always nine electrons in 
the impurity d orbitals. The temperature is k^T — 0.01 
eV as before. 
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FIG. 5. (color online) Spin-resolved d-orbital spectral function 
of a cluster with parameters taken from Ref. and tb = 
eV (lower panel). The "envelope" is calculated with a large 
Lorentz broadening of 0.5 eV. The top panel shows the d- 
orbital occupation in the final states. 

The spin polarization of the impurity orbitals is induced 
by a small magnetic field B coupled to the impurity spins. 
The coupling is described by an extra term in the cluster 
Hamiltonian, 

Hi""^ - f E (^lt^™t - dl^d„^i) ■ (Bl) 

m 

The resulting polarization is plotted in Fig.[4]as a function 
of the bath position eb scanned across the Fermi level. 
The total number of electrons in the cluster is shown as 
well. 

When the bath orbitals are sufficiently deep below the 
Fermi level, the bath is nearly full and a local magnetic 
moment is formed on the impurity. As the bath orbitals 
move up toward the Fermi level, the bath relatively quickly 
depopulates until it holds only a single electron. This 
electron, together with the other nine sitting in the im- 
purity orbitals, forms a non-magnetic d}^ closed shell. 
This state then remains stable even when the bath is 
raised relatively high above the Fermi level. Only for 
Eb > 2.4VfeTOcr ~ 1-7 eV the bath starts releasing the last 
electron and a magnetic ground state is restored. The 
larger cluster corresponding to our actual DMFT calcula- 
tions shows an analogous behavior, only the non-magnetic 
solution occurs for 20 electrons in the cluster as there 
is an extra fully occupied shell of bath orbitals located 
deeper below the Fermi level. 

It is clear that the non-magnetic solution does not 
correctly describe the d shell and its environment in the 
ferromagnetic nickel. Elevating the bath orbitals high 
above the Fermi level in order to support a magnetic 
ground state does not look as very plausible. This leaves 
us with the configuration where the bath states are nearly 
fully occupied and hence they model the nearly full d 
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orbitals of the nickel atoms surrounding the "impurity" 
site. To prevent the cluster Hamiltonian to enter the non- 
realistic regimes in the course of the DMFT iterations, 
we have introduced the upper cutoff in the integrals in 
Eq. ([9| . This cutoff does not allow the bath orbitals to 
drift too high and to lock into the non-magnetic solution. 

It is instructive to compare the spectral functions corre- 
sponding to the different cluster ground states. Figure [5] 
shows the spectrum obtained when the bath orbitals are 
placed exactly at the Fermi level, Cb = 0. The local mo- 
ment induced by the external magnetic field is negligible 
in this case and the spectral function is nearly symmet- 
ric. The spectrum is practically identical to the result 
presented in Ref. JJj as it should be, since we used the 
same parameters. The superimposed plot of the d-orbital 
filling nj^ in the photoemission final states indicates that 
the satellite structures around —6.5 eV and near —10 eV 
are due to the final states. 

The spectral function corresponding to the bath orbitals 
lowered to —0.4 eV is plotted in Fig.|6| The "6 eV satellite" 
has now a shape similar to our DMFT solution (Fig. ^ 
as well as to the experimental data: the minority-spin 
component is less intense and is located at smaller binding 
energies. Comparison of Figs, [s] and [6] reveals that the 
composition of the final states constituting the main d 
band is shifted toward a larger average n^, likely due to 
an increased impurity-bath hybridization caused by the 



reduced distance between the impurity and bath orbitals. 
Analogously, the enhanced values Ud ~ 8.5 at the satellite 
in the DMFT solution (Fig. [2]) are probably due to the 
hybridization with the extra shell of bath orbitals not 
present in the model discussed in this appendix. 
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FIG. 6. (color online) Spectrum of the same model as Fig. [s] 
but with the bath orbitals shifted to —0.4 eV. The ground 
state is magnetic in this case. Note the different range of the 
top panel compared to Fig. [5] 
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